Observation of Three-dimensional Long-range Order 
in Smaller Ion Coulomb Crystals in an rf Trap 
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Three-dimensional long-range ordered structures in smaller and near- spherically symmetric 
Coulomb crystals of 40 Ca + ions confined in a linear rf Paul trap have been observed when the num- 
ber of ions exceeds ^1000 ions. This result is unexpected from ground state molecular dynamics 
(MD) simulations, but found to be in agreement with MD simulations of met ast able ion config- 
urations. Previously, three-dimensional long-range ordered structures have only been reported in 
Penning traps in systems of ^50,000 ions or more. 

PACS numbers: 32.80.Pj, 52.27.Jt, 52.27.Gr, 36.40.Ei 
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A Coulomb crystal is the solid state phase of a con- 
fined ensemble of Coulomb interacting particles with the 
same sign of charge, often referred to as a one-component 
plasma (OCP). Coulomb crystallization and properties 
of Coulomb crystals have been studied experimentally 
for decades in such a variety of systems as 2D electron 
gases on super- fluid helium [l| and in quantum well struc- 
tures , laser cooled ions in traps 0, El IS El El 
[H EH and most recently dusty plasmas [TJ. In nature, 
Coulomb crystals are presently expected to exist in exotic 
dense astrophysical objects 

Theoretically, it has been found that the thermody- 
namic properties of infinite OCPs of a single species, are 
fully characterized by the coupling parameter |T| 
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where Q is the charge of the particles, a ws is the Wigner 
Seitz radius denned from the zero temperature particle 
density no by 47ra^ s /3 = I/uq. Furthermore, a liquid- 
solid transition to a body centered cubic (bcc) structure 
is expected to occur for T ~ 170 0, For finite 

OCPs simulations have shown that the situation is more 
complex, and the properties will depend both on the size 
and shape of the ion plasma , si nce surface effects cannot 
be neglected [H El 111 E3 

Ion Coulomb crystals, which for more than a decade 
have been realized with laser-cooled ion plasmas confined 
by electromagnetic fields in Penning trap s 0, H or in 
radio- frequency (rf) traps 0, H 0, H llOj . offer an excel- 
lent opportunity to study finite size effects of OCPs un- 
der various conditions. The Coulomb crystal structures 
studied range from one-dimensional (ID) long cylindrical 
crystals 0000 over 2D thin planar crystals to 3D 
spheroidal crystals 0, 0, H 0. The 3D spheroidal ion 
Coulomb crystals reported in Refs. are composed 
of concentric ion shells formed under the influence of the 
surface of the Coulomb crystals. Simulations indicate 
that the ions form a near-2D he xag onal short-range or- 
dered structure within each shell PJ| . Similarly, shell and 



short-range order have recently been observed in dusty 
plasma experiments [13 . 

Observations of three-dimensional long-range order in 
Coulomb crystals have previously only been reported in 
the case of ^50,000 or more laser-cooled ions in a Pen- 
ning trap 000. In contrast to Penning traps, in rf 
traps Coulomb crystals undergo strong quadrupole de- 
formations at the frequency of the applied rf field due to 
the so-called micro- motion |22j of the ions. Since this mo- 
tion is known to produce heating [21 [2^ , it has not been 
obvious that three-dimensional long-range order could be 
obtained in such traps. 

In this Letter, we present observations of long-range 
structure in Coulomb crystals of 40 Ca + ions confined in 
a linear rf Paul trap. By varying the number of ions in 
near-spherically symmetric crystals, we have shown that 
bcc structures indeed can be observed in such traps even 
with the number of ions being below a thousand. 

The linear Paul trap used in the experiments is de- 
scribed in detail in Ref. 24], so here we give only a 
brief description of its basic properties. The linear Paul 
trap consists essentially of four circular electrode rods in 
a quadrupole configuration. Radial confinement of the 
ions is obtained by applying a sinusoidally time-varying 
rf-potential U T f cos Q r ft to diagonally opposite electrodes 
and the same time varying potential, but with a phase 
shift of 7r, to the two remaining electrode rods. Static ax- 
ial confinement is achieved by having all electrode rods 
sectioned into three parts with the eight end-pieces kept 
at a positive dc potential U en d with respect to the four 
center pieces. The electrode rod radius is 2.0 mm, and 
the closest distance from the trap axis to the electrodes ro 
is 3.5 mm. Each sectioned rod has a 5.4 mm long center- 
piece and 20.0 mm long end-pieces. The frequency of the 
rf field is in all experiments O^r = 2tt x 3.88 MHz. With 
the definition given in Ref. |23| of the relevant stability 
parameters q and a for the linear Paul trap under the 
present operation conditions, the geometry of the trap 
leads to q = 6.6 x 10- 4 Z7 rf [V] and a = 5.5 x 10- 4 U end [V]. 
In the experiments the rf-amplitude U Y { [V] is maximally 
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500 V, corresponding to q max = 0.33. This value is small 
enough that in all the experiments a harmonic pseudopo- 
tential with rotational symmetry with respect to the trap 
axis and a uniform zero temperature ion density given by 
ntheo = e /7 r 2 f /rar^ 2 f = L47 x io 3 /7 rf [V] 2 cm" 3 can be 
assumed (eo is the vacuum permittivity and m the mass 
of the ion). 

The 40 Ca + ions are produced isotope-selectively by 
resonance-enhanced photo-ionization of atoms from an 
effusive beam of naturally abundant calcium H3- 
Doppler laser cooling is achieved along the trap center 
axis by using counter propagating laser beams tuned to 
the 46 , 1 / 2 -4P 1 /2 transition at 397 nm and with a re- 
pumper beam on the 3D 3 / 2 ~^Pi/2 transition at 866 nm. 
By this choice of cooling geometry the fluorescence level 
of the individual ion will not significantly be effected by 
the rf-induced micro motion (see Fig. 1), as this motion 
is perpendicular to the direction of the cooling laser to 
better than 1 mrad. The radial motion of the ions is 
cooled indirectly through the Coulomb coupling of the 
radial and axial degrees of freedom. 

Images of the positions of the ions in the trap are ob- 
tained by detecting the 397 nm light spontaneously emit- 
ted during the laser cooling process by a CCD-camera 
equipped with an image intensifler and a 14 x magnifica- 
tion lens system. As indicated in Fig. 1(a), the imaging 
optic is situated such that the real crystal structure is 
projected to a plane including the trap axis. Due to the 
rotational symmetry of the Coulomb crystal boundary 
around the trap axis, we can deduce from the images 
their real 3D sizes and hence the number of ions by using 
the expression for the ion density given above. 

With an exposure time of 100 ms used in the experi- 
ments, the quadrupole deformations of the Coulomb crys- 
tal induced by the rf field are averaged out in the images. 
The direction of the micro-motion of the individual ions 
is position dependent as indicated in Fig. 1(a). As a re- 
sult, only ions close to the horizontal trap plane defined 
by the horizontal dashed line in Fig. 1(a) and the trap 
axis, are imaged without micro-motion blurring due to 
their micro-motion only being in the direction of view 
(DV) of the camera system. 

In Fig. 1(b), an image of a slightly prolate ion Coulomb 
crystal consisting of about 2300 ions is presented. Clearly 
visible is a hexagonal structure. Since the depth of fo- 
cus (~50 /im) is several times the inter-ion distance (~15 
/im), the observed structure cannot originate from a sin- 
gle layer of ions. In order to reveal if the ordered struc- 
ture persists not only very close to the image plane, 
in another experiment, the Coulomb crystal was moved 
down a few inter-ion distances (~40 /im) by adding a 
positive dc voltage to the upper two electrode rods as 
sketched in Fig. 1(c). By doing this one achieves that 
it is a different part of the crystal that is close to the 
horizontal trap plane. As seen in the image presented 
in Fig. 1(d), in this situation we still observe an ordered 




FIG. 1: Observation of hexagonal structures in images of a 
Coulomb crystal consisting of ^2300 ions, (a) The four main 
circles indicate the four cylindrical trap electrode rods viewed 
along the trap axis. The small circle in the center represents 
a radially centered Coulomb crystal in the trap. The small 
double arrows indicate the positional dependent direction of 
the micro-motion of the ions in the trap. The single arrow 
denotes the direction of view (DV) for the camera system. 
The horizontal dashed line indicates the plane including the 
trap axis in which the direction of micro-motion is along DV. 
(b) Image with the Coulomb crystal radially centered in the 
trap, (c) Graphical presentation of a crystal vertically dis- 
placed, (d) Image of a vertically displaced crystal. In both 
experiments: [7 r f=400 V and £/ en d=15 V. 



structure apart from close to the edges where the flu- 
orescence light originates from ions not being close to 
the image plane. This experiment substantiates that the 
ordered pattern observed indeed originates from a three- 
dimensional long-range ordered structure. Furthermore 
from a density point of view the observed structure can 
not just be two-dimensional. 

From the pre-knowledge that the grou nd state of larger 
ion systems is a bcc configuration jlHl H(| , it is tempt- 
ing to assume that the observed structure is indeed a 
bcc structure observed along the [111] direction. How- 
ever, also simple cubic (sc) and face centered cubic (fee) 
structures will lead to similar hexagonal patterns for the 
same direction of observation. For a specific side length 
d of the triangles making up the observed structure, the 
three cubic structures correspond, however, to different 
ion densities related by T^bcc = 

2n sc = 4n fcc = 1.089/d 3 . 
The perfect agreement between the expected density of 
^theo = 2.3 ± 0.2 x 10 8 cm -3 with the one obtained from 
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the images nb cc = 2.1 ± 0.3 x 10 8 cm -3 assuming a bcc 
structure, strongly supports that the observed structure 
is indeed a bcc and certainly not a sc or fee structure. 
As seen from Fig. 2(a)-(c), the outer shape of the crys- 
tal is not critical for observing the long-range structures, 
at least not as long as the crystal is not too oblate or 
prolate. 




FIG. 2: Images of Coulomb crystals. The first column shows 
three crystal images of the same ensemble of ^2000 ion with 
the trap potentials being U T f = 500 V and (a) U en d = 27 V, 
(b) Uend = 33 V, and (c) U en d = 40 V, respectively. The 
second column represents images near-spherical crystals for 
the same trap potentials U T f — 400 V and U en .d = 21 V. The 
number of ions in the crystals is (d) 1700, (e) 770, and (f) 
290. 

A series of measurements on near spheric crystals with 
different numbers of ions were done to investigate the 
size effect of the observation of three-dimensional long- 
range structures. A few resulting images of this study 
are presented in Fig. 2(d)-(f). Here one observes long- 
range-order for as few as 770 ions (Fig. 2(e)), and even 
some structure is present in the core for the case of 290 
ions (Fig. 2(f)). 

When comparing the images in Fig. 2, it is evident 
that the relative number of ions in the regular struc- 
ture decreases with the size of the crystals, as may be 
expected since the outer layer of the crystals is always 
spheroidal shaped. Another size dependent quantity no- 
ticed is the frequency at which ordered structures are 



indeed observed: The smaller the crystals the more un- 
likely it is to observe the regular structures. The rate at 
which the long-range ordered structures are observed in 
crystals of a few thousand ions is ~0.1 Hz and the life- 
time of the structure is a few hundred ms. In all cases, 
since the exposure time of 100 ms is many orders of mag- 
nitude larger that both the rf-period and the timescale 
of crystal vibrations (~l/c^piasma? with cj p i a sma ~ 

1 MHz 

being the plasma frequency), the crystal structures are 
at least to be considered as metastable states. 

Molecular dynamics (MD) simulations of ground state 
configurations of Coulomb crystals in static spherical 
harmonic potentials have proven only to have a pro- 
nounced bcc structure when the number of ions exceeds 
~5000 [^ . E^ . However, in contrast to such ground state 
simulations, a certain thermal energy is always present 
in experiments, and configurations different from the 
ground state may occur. In order to understand the ob- 
served long-range-ordered structures, we have made sev- 
eral MD simulations using a static spherical harmonic po- 
tential where we initially kept a core of the ions fixed in a 
bec-structure and cooled the surrounding ions slowly un- 
til T > 100.000 was reached. For an ion system of ~1000 
ions, the excess potential energy (relative to the ground 
state configuration) of such artificially created cold con- 
figurations is even with more than one tenths of the ions 
kept fixed in a bec-structure, several times smaller than 
the thermal energy of the system at the typical experi- 
mental temperature of a few mK. Simulations where all 
the ions were subsequently slowly heated to F ~ 400 (a 
few mK) furthermore showed that the artificially made 
configurations could be metastable on the time scale of 
the exposure time of the images, ~10 ms. Figure 3(a) is 
an example of a constructed image based on the results 
of a MD simulation of a crystal of 2685 ions heated to 
T = 400 (temperature: 5 mK), where initially 10% of 
the ions were kept fixed in a bec-structure. With the 10 
ms integration time of this simulation, the central bec- 
structure is seen to be well preserved, and the image re- 
sembles very much the experimentally obtained image of 
a Coulomb crystal with about the same number of ions 
(Fig. 3(b)). 

For crystals with more than ^2000 ions, images are 
observed, which suggests three-dimensional long-range 
ordering different from bcc. In Fig. 4, two images of 
the same ion ensemble (~13,000 ions) at two different 
instants are presented. While the hexagonal structure 
of Fig. 4(a) and (b) with d = 17.2 ± 0.8 /im again is 
compatible with a bcc structure, the rectangular struc- 
ture observed in Fig. 4(c) and (d) cannot be interpreted 
as another projection of a bcc structure, but has to re- 
late to another structure. The sides of the rectangle are 
h = 15.2 ± 0.5 /im and w = 9.7 ± 0.2 /im, which yields a 
ratio of the sides of h/w = 1.55 ± 0.06. We do not find 
that any projection of cubic crystal structures exactly 
complies with this ratio, but the fee structure projected 




FIG. 3: Images of Coulomb crystals, (a) Time averaged im- 
age based on data from MD simulations of Coulomb clusters 
with 2685 ions at F ~ 400 (temperature: ^5 mK). The aver- 
aging time is 10 ms. (b) Image from experiments with clusters 
containing ~2700 ions. 




FIG. 4: Images of a cold ensemble of ~ 13,000 ions, (a) Visible 
hexagonal structures indicating a three-dimensional bcc struc- 
ture, (b) Magnification of a section of (a). The side length 
d is found to be 17.2 /xm. (c) Visible rectangular structure 
likely indicating a slightly distorted fee structure, (d) Mag- 
nification of a section of (c). Here, h = 15.2 db 0.5 /mi and 
w = 9.71 ±0.2 jLtm, respectively. 



in the [211] direction comes close as this has the ratio 
h/w = \/S/3 ~ 1.63. For the expected ion density of 
^theo = (2.2 ± 0.2) x 10 8 cm -3 the corresponding side 
lengths of the projected rectangles should for a fee struc- 
ture have been h = 15.2 ± 0.5 /im and w = 9.3 ± 0.3 jam. 
The origin of the small discrepancy between the observed 
and the predicted ratio is at present unknown, but an ex- 
planation could probably be found in a micro-motion in- 
duced distortion of the fee lattice. The observation of fee 
structures is not so surprising since such structure have 
already been observed for larger near-spherical Coulomb 



crystals in Penning traps 5]. Furthermore, MD simu- 
lations have as well predicted very small differences in 
the potential energies of bcc and fee structures for larger 
crystals |29j . 

The reason why the observed bcc or fee structures al- 
ways seem to observed from a specific direction is still 
not completely understood. However, the non-perfect 
cylindrical symmetry of the quadrupole trap configura- 
tion and the corresponding micro-motion as well as small 
patch potentials may play a role. 

In conclusion, we have proven that it is possible to 
obtain three-dimensional long-range ordering in smaller 
ion Coulomb crystals formed in rf traps. The long-range 
structures are found to be consistent with metastable 
structures observed in MD simulations. 
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Foundation: Center for Quantum Optics QUANTOP. 
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